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Abstract 

The regions in the Minimal Supersymmetric Standard Model 
with the minimal amount of fine-tuning of electroweak symmetry 
breaking are presented for general messenger scale. No a priori 
relations among the soft supersymmetry breaking parameters are 
assumed and fine-tuning is minimized with respect to all the im- 
portant parameters which affect electroweak symmetry breaking. 
The superpartner spectra in the minimally tuned region of pa- 
rameter space are quite distinctive with large stop mixing at the 
low scale and negative squark soft masses at the high scale. The 
minimal amount of tuning increases enormously for a Higgs mass 
beyond roughly 120 GeV. 
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1 Introduction 

The Minimal Supersymmetric Standard Model (MSSM) is a well-motivated 
candidate for physics beyond the Standard Model (SM). The gauge couplings 
within the MSSM unify to within a few percent at the grand unified theory 
(GUT) scale, M GUT ~ 2 x 10 16 GeV, and the lightest supersymmetric par- 
ticle is a good dark matter candidate provided that R-parity is conserved. 
Supersymmetry (SUSY) can also naturally stabilize the hierarchy between 
the electroweak (EW) and the GUT or Planck scale. It does this by pro- 
viding a radiative mechanism for electroweak symmetry breaking (EWSB) 
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where large quantum fluctuations of the scalar top squarks due to the large 
Yukawa coupling destabilize the origin of the Higgs potential. In much of the 
MSSM parameter space this quite naturally leads to the right EWSB scale, 
as long as the soft SUSY breaking parameters lie near it. 

The absence of any direct experimental evidence from collider searches 
for the MSSM scalar particles and the Higgs boson has, however, ruled out 
significant regions in the MSSM parameter space. Indirect evidence from 
EW precision measurements and searches for flavor changing neutral cur- 
rents, CP violating effects and rare decays has not been forthcoming either, 
providing additional severe constraints. As a result the soft SUSY breaking 
parameters must lie well above the EW scale in order to satisfy the exper- 
imental constraints, especially the constraints on the Higgs mass from the 
results of the CERN LEP collider (m h > 114.4 GeV [1]). 

Soft SUSY breaking parameters well above the EW scale reintroduce a 
small hierarchy and require some fine-tuning (FT) among the SUSY param- 
eters in order to obtain EWSB [2j- [22]. This is usually referred to as the 
supersymmetric little hierarchy problem. 

Different choices for the soft SUSY breaking parameters lead to differ- 
ent amounts of FT. This paper presents the minimally tuned MSSM (or 
MTMSSM), i.e. the MSSM parameter region that has the least model- 
independent FT of EWSB. Model- independent means that no relations are 
assumed between the soft SUSY breaking parameters at the scale at which 
they are generated (which will be referred to as the messenger scale) . Rather, 
each of them is taken to be an independent parameter which is free at the 
messenger scale, and which therefore can contribute to the total FT of the 
EWSB scale. The messenger scale itself is varied between 2 TeV and Mqut 
and the effect of this on the minimal FT is discussed (see also [141 [T5]). 

In Section [2J, EWSB in the MSSM will be reviewed. Section [3] discusses 
the tuning measure used in this paper. The parameters taken to contribute 
to the tuning are |/i| 2 , m 2 Hu , the gaugino masses Mi, M 2 and M 3 , the stop 
soft masses m~ and m~ , and the stop soft trilinear coupling A t . 

Section [4] contains some of the main results. The low- and high-scale 
MSSM spectrum which leads to the least model-independent FT is found. 
This is done for various messenger scales by numerically minimizing the FT 
expression subject to constraints on the Higgs, stop, and gaugino masses. The 
results are then motivated analytically. The least FT is found to be about 5% 
if the messenger scale coincides with the GUT scale. An important feature 
of the least FT region is negative stop soft masses at the messenger scale 
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(first pointed out in [IS])- Even for messenger scales as low as 2 TeV, the 
stop soft masses are tachyonic at the messenger scale (threshold effects in the 
RG-running were neglected throughout). This does not lead to any problems 
with charge and/or color breaking minima. Another feature of the least FT 
region is that the trilinear stop soft coupling, A t , is negative and lies near 
"natural" maximal mixing, i.e. A t ~ — 2m t ~, where is the average of the 
two stop soft masses. This value for A t maximizes the radiative corrections to 
m/j. The large stop mixing leads to a sizeable splitting between the two stop 
mass eigenstates. Moreover, the gluino mass, M3, is much smaller than the 
wino mass, M2, at the high scale. The wino mass, in turn, is much smaller 
than the bino mass M\. Phenomenological consequences of the low-scale 
spectrum are briefly summarized. 

Section [5] contains the rest of the main results of the paper. The FT 
is minimized as a function of the lower bound on the Higgs mass (with 
the messenger scale set to Mgut)- Although the numerical minimization 
procedure contains the dominant one- loop expression for as a constraint, 
the resulting least FT spectra are used to calculate m h more accurately with 
the program FeynHiggs [231 EU E51 ESI EH- The result is a plot of the 
minimal FT as a function of mh- There are several striking features of this 
plot. First of all, for larger than a certain value, the FT increases very 
rapidly. Secondly, around this m^, the value of A t in the least FT region 
makes a sudden transition from lying near — 2m^ to lying near +2m^. The 
third striking feature is that this value of is surprisingly low. The precise 
value is only slightly dependent on the parameters in the Higgs sector and can 
be taken to lie around 120 GeV. The upshot of this analysis is that although 
the MSSM right now is already fine-tuned at least at about the 5% level (if 
the messenger scale equals the GUT scale), there is not much room left for 
the Higgs mass to increase before the FT becomes exponentially worse. 

Section [6] contains a summary of the results and the conclusions. Ap- 
pendix [A] reviews the semi-numerical solutions of the MSSM one-loop renor- 
malization - group (RG) equations. These are used to calculate the expression 
for the FT employed in this paper. Appendix [B] contains a list of expressions 
for the FT with respect to various parameters. 
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2 Electroweak Symmetry Breaking 



In the Higgs decoupling limit of the MSSM, the lower bound on the mass of 
the lighter CP-even Higgs mass eigenstate h coincides with the 114.4 GeV 
bound on the mass of the SM Higgs boson pQ. The mass of h may be 
approximated by 



m 2 h ~ m 2 z cos 2 2(3 + — -f 



m i , / X 2 



log -4 + -4 



m 2 m~ \ 12m 2 



which, in addition to the tree-level Higgs mass, includes the dominant one- 
loop quantum corrections coming from top and stop loops EDI ED E2J 
S3] . Here m t is the top mass, m~ is the arithmetic mean of the two squared 
stop masses and v = 2mw/g — 174.1 GeV where g is the SU{2) gauge 
coupling and my/ is the mass of the H^-boson. Furthermore, equation (pP) 
assumes m~ t ^> m t . The stop mixing parameter is given by X t = A t — fi cot (3 
(~ At for large tan/3), where A t denotes the stop soft trilinear coupling 
and n is the supersymmetric Higgsino mass parameter. The first term in 
equation (CTJ is the tree-level contribution to the Higgs mass. The first term 
in square brackets comes from renormalization group running of the Higgs 
quartic coupling below the stop mass scale and vanishes in the limit of exact 
supersymmetry. It grows logarithmically with the stop mass. The second 
term in square brackets is only present for non-zero stop mixing and comes 
from a finite threshold correction to the Higgs quartic coupling at the stop 
mass scale. It is independent of the stop mass for fixed X t /mi, and grows as 
(X t /mi) 2 for small X t /mf. 

Equation ([1]) implies a combination of three things which are required to 
satisfy the bound on m^, namely a large tree-level contribution, large stop 
masses and large stop mixing. A large tree-level contribution to m^ requires 
tan/3 to be at least of a moderate size (> 5 — 10). Although the stop masses 
must be rather large, their lower bound is very sensitive to the size of the 
stop mixing, with larger mixing allowing for much smaller stop masses (see 
[33] for a recent study on this). The reason for this sensitive dependence 
is due to the Higgs mass depending logarithmically on the stop masses in 
contrast to the polynomial dependence on the stop mixing. 

The soft masses are not only directly constrained from the LEP Higgs 
bounds but also indirectly by constraints on flavor changing neutral currents, 
electroweak precision measurements and CP-violation. Besides these, how- 
ever, the Higgs sector parameters are also constrained by requiring that the 
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electroweak symmetry is broken. This leads to the following two tree-level 
relations at the low scale 

sin 2/3 = 2 2 f 2 . = (2) 



TJlr 



mir — mir tan 8 

rid tlv. " 



2 1/11 + tan 2 /3-l ' [6) 
where tua is the CP-odd Higgs mass, and 8 is determined from the ratio 
of the two vacuum expectation values v u = (Re(if°)) and Vd = (Re(if°)) 
as tan/3 = v u /v<i- The masses tt^, Jnfr, and are the three soft mass 
parameters in the MSSM Higgs sector. For a given value of tan 8, m 2 2 may be 
eliminated in favor of m 2 A with equation ([2]). Equation ([3]) gives an expression 
for m| in terms of the supersymmetric mass parameter /i and the soft masses 
m 2 H and m 2 H . Since tan 8 should be sizeable, the contribution from m 2 H to 
the expression for m 2 z may be neglected and ([3]) simplifies to 

m| = -2|/i| 2 -2m^. (4) 

Close to the Higgs decoupling limit, rriA is relatively large. However, since 
\/j,\ 2 ,m 2 Hu ~ 0(m 2 z ) to avoid large cancellations, vtla may not be too large, 
otherwise rn 2 H would also be sizeable and equation (j3J) would break down 
(unless the value of tan/3 is increased accordingly). By choosing tan/3 = 10 
and rriA = 250 GeV in the numerical analysis throughout, equation (TJJ holds 
to a very good approximation. 

Equation (TjO) holds at tree-level, and although quantum corrections may 
add O(10 GeV) to the right hand side of (jl]), this has negligible impact on 
the amount of fine-tuning to be discussed below. 

The parameters rn 2 H and |/z| 2 in equation (j4j) are evaluated at the scale 
mz- Since the fine-tuning of EWSB is a measure of the sensitivity of some 
low-scale EWSB parameter (usually taken to be m 2 z ) to a change in high- 
scale input parameters, |/i| 2 and m 2 H need to be evolved to a high scale using 
their RG equations. Under RG running many of the soft parameters mix, and 
as a result of this mixing, the expression for m 2 z in terms of parameters that 
are evaluated at the messenger scale Ms differs significantly from the simple 
form given in (jl]). The RG-equations may be integrated (see Appendix 1X1) 
and the expression for m 2 z may generically be written as 



^ % (tan/3, M s ) m^M^m^M^). (5) 
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Fig. 1: The coefficients Cij defined in equation © for tan/3 = 10 as a function of 
the messenger scale Mg. 



For moderate and not too large values of tan/5 with an appropriate m^, the 
simplified expression for ra| is applicable (equation and contributions 
from the bottom/sbottom and tau/stau sectors may still be neglected. The 
most important parameters appearing in (jSJ) then are /i 2 , m 2 H , the gaugino 
masses Mi, M2 and M3, the stop soft masses m~ and m~ , and the stop soft 
trilinear coupling A t . The coefficients depend on tan/5 and the messenger 
scale Ms- The most important coefficients are shown in Figured] for tan/3 = 
10 as a function of Mg. 

At the scale mz, the coefficients of mf^ and /x 2 are —2 while the coeffi- 
cients of the other soft parameters are zero in agreement with equation (j3J). 
Since fi 2 is a supersymmetric parameter, it gets renormalized multiplicatively 
and its RG evolution does not give rise to soft parameters (see equation f|45|) ). 
Figure [1] shows that the coefficient of /x 2 does not vary much and remains 
close to —2 all the way up to the GUT scale. The RG evolution of rn 2 Hu 
to higher messenger scales, however, generates non-zero coefficients for the 
other soft parameters. The /3-function of m 2 H , 

87r 2 /3 m 2 Hu =3A t 2 (m^+m^^^ (6) 
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depends on the stop sector parameters {m~^,m 2 ^, A t }, the wino and bino 
masses M 2 and Ml, and Sy = |Tr(y^rnf), which thus get generated im- 
mediately under RG evolution. The coefficients of M 2 and especially Mi 
and Sy in ([6]) are small and lead to small coefficients in the expression for 
m| flSJ). Although f3 m 2 does not explicitly depend on the gluino mass, a 
non-zero coefficient for M3 is generated indirectly since the stop sector (3- 
functions depend on M3. Moreover, M3 appears with a large coefficient in 
these /3-functions, and thus the coefficient of M3 in equation (jSJ) dominates 
after a few decades of RG evolution. For example, at a messenger scale of 
Ms = Mqut = 2 x 10 16 GeV, the expression for m| (for tan/3 = 10) is 

m| = —2.19 /i 2 - 1.32 m 2 ^ + 0.68 m?. + 0.68 m 2 ^ + 5.24 M 3 2 - 0.44 M 2 
-0.01M 2 + 0.22 A 2 - 0.77 A t M 3 - 0.17 A t M 2 - 0.02 i t M x 
+ 0.46 M 3 M 2 + 0.07 M 3 Mi + 0.01 M 2 Mi + 0.05 Sy, (7) 

where the hatted parameters on the right-hand side are all evaluated at Ms. 
This expression may be used to calculate the FT as discussed next. 



3 The Tuning Measure 



A variety of tuning measures have been used in the literature (a list of refer- 
ences has been provided in the Introduction). Since the concept of fine-tuning 
(FT) is inherently subjective, there is no absolute definition of a FT measure. 
The most common definition of the sensitivity of an observable 0{{a{\) on 
a parameter denoted by A(C, a^), is given by [2l [3] 



A(0, ai 



dhgO 




ai dO 


d log at 




Oda~i 



(8) 



A((9,aj) thus measures the percentage variation of the observable under a 
percentage variation of the parameter. A large value of A(0, a;) signifies that 
a small change in the parameter leads to a large change in the observable, 
and suggests that the observable is fine-tuned with respect to that param- 
eter. Assuming that the individual A(0,ai) are uncorrelated, they may be 
combined to form the FT measure 



T{0) 




(9) 



Of interest in this paper is to quantify the sensitivity of EWSB in the 
MSSM on (soft) supersymmetric parameters at the messenger scale Ms- To 
this end, the observable to consider is m| as a function of the supersymmetric 
Higgsino mass squared and the soft supersymmetry breaking parameters, 
collectively denoted by mf(Ms) (in the FT measure, all parameters are taken 
to have mass dimension two). The sensitivity of m| with respect to each 
parameter may be calculated as in ([S]) with O = m|, and the total FT of m| 
on parameters evaluated at the messenger scale Ms may be quantified by 

H<\ Ms) = (A(m|X(M 5 ))) 2 . (10) 

jF(m|; Ms) may be interpreted as the length of a "fine-tuning vector" with 
components A(m%, mf (Ms))- This fine-tuning vector is formally a vector 
field defined by the gradient of the scalar field logm|, a function of logmf, 
along surfaces of constant logmf-. 

There are several possible drawbacks to this FT measure, see for example 
[22l I3T] . One of these is that the individual A(m|, mf (Ms)) are assumed 
to be uncorrelated. Within a given model of supersymmetry breaking, there 
may be relations among the parameters at the messenger scale. This would 
imply that the FT vector is projected onto a subspace, and the resulting FT 
is necessarily less. In other words, the tuning of one parameter is correlated 
with the tuning of another, so that the total FT should be less than that given 
by ffTUj) . Moreover, within a given model the values of the parameters at the 
messenger scale may be restricted to certain ranges, whereas f fTUj) assumes 
that all values are equally likely. However, no model for supersymmetry 
breaking will be assumed here. Instead, the minimal FT will be found as 
a function of the messenger scale Ms assuming no relations or restrictions 
among the high-scale input parameters. For this "model-independent" tuning 
it is satisfactory to use the FT measure (fTUj) . 

Note that to find the tuning of a model, one should in principle consider 
the tuning of all observables, since the absence of tuning in one observable 
does not necessarily imply it is small in others, see e.g. [T7]. In this paper, 
however, only the tuning of EWSB will be considered. 

Finally, note that the FT with respect to a single parameter is by defini- 
tion (J8]) zero if that parameter happens to be zero at the messenger scale. An 
extreme version of this is found in the no-scale model [38], where all scalar 
soft masses are much smaller than the gaugino masses at the high scale. 
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Setting them to zero, and using (jSJ) and (TTOT) the FT could be expected to 
be small. However, it may be shown that this does not minimize the FT, 
since M 3 and /i need to be quite large at the high scale to satisfy all the 
low-energy experimental bounds (see [13]). In the results presented in this 
paper, no parameter is found to be zero at the high scale. 

4 Minimal Model Independent Tuning 

In this section the minimal model independent tuning will be found as a 
function of the messenger scale. 

4.1 Discussion of Minimization Procedure and Con- 
straints 

The FT given by equation fflOl) is written in terms of parameters evaluated 
at the messenger scale. In order to find the minimal FT (MFT) for a given 
messenger scale that is consistent with low-energy experimental constraints, 
it is easiest to rewrite the FT expression in terms of parameters that are 
evaluated at the low scale. This can be done by expressing each high-scale 
parameter in terms of low-scale parameters, see Appendix |A] Once the FT 
is written in terms of low-scale parameters, m 2 Hu (mz) may be eliminated by 
using equation (J4]) (neglecting contributions from m 2 H ). 

The low-energy constraints considered in this paper include bounds on the 
(physical) sparticle masses, on the gaugino masses, and on the Higgs masfl 
The top quark mass m t is set to the central value of the latest Tevatron mass 
measurement of 170.9 ± 1.8 GeV [39]. The physical stop masses are required 
to be at least 100 GeV which is illustrative of the actual, slightly model 
dependent, lower bound obtained from the Tevatron [10]. It is found that 
the region of MFT does not quite saturate this bound, although a slightly 
larger value for m t would allow the lighter stop to be as low as 100 GeV. The 
gaugino masses M\ and M2, as well as fi, are taken to have a lower bound 
of 100 GeV. The gluino mass is found to be never smaller than 335 GeV in 

1 Constraints from measurements of B — > X s j or the electroweak S- and T-parameter 
do not significantly affect the results presented below, since an experimentally consistent 
value can be obtained by only small adjustments (if at all necessary) in the least fine-tuned 
parameters - see also [34] . 
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the numerical results presented in this section, and this does not generically 
violate any experimental bounds. 

The most important constraint is the Higgs mass bound of 114.4 GeV 
(valid in the decoupling limit), since it turns out that this bound is always 
saturated when minimizing the FT. In the numerical results presented in 
this paper, the Higgs mass is calculated using the formulas found in [JT| (see 
also [221 EH ED E21 E2] ) • These formulas include the one-loop corrections 
coming from the top/stop sector and are simple enough to be used as con- 
straints in the FT minimization (but note that the sign convention used for 
A t is that of [13]). A running top mass m t (m t ), evaluated in the MS'-scheme, 
is used to capture some of the leading two-loop contributions. Higher-order 
corrections to the Higgs mass still play a very important role, however, and 
more accurate Higgs masses may be obtained with the program FeynHiggs. 
In order to take some of these higher-order corrections into account and thus 
obtain a more accurate estimate of the MFT, a lower bound for the Higgs 
mass of 121.5 GeV is used in the FT minimization, instead of the SM lower 
bound of 114.4 GeV. The typical low energy sparticle spectrum obtained in 
the analysis then leads to Higgs masses that satisfy the SM Higgs bound 
when calculated with FeynHiggs (version 2.6.0, assuming real parameters). 

Sequential Quadratic Programming (SQP) is used as a minimization algo- 
rithm. Given the FT function (fit)]) written in terms of low scale parameters, 
as well as linear constraints on the gaugino masses and fi, non-linear con- 
straints on the physical stop and Higgs masses, and an initial guess, SQP 
generates a less FT point until the minimum is found. Unlike other mini- 
mization algorithms, SQP can handle arbitrary constraints which is essential 
here due to the highly non-linear physical stop mass and Higgs mass con- 
straints. 

4.2 Numerical Results 

Figure [2] shows a plot of the MFT as a function of the messenger scale Ms- 
Shown are the individual contributions A(m|, mf (M^)) to the FT, with m] 
given by Mf , Mf, Mf, A 2 , /x 2 , or m 2 Hu . The FT of m 2 ^ and m 2 ^ have been 
included as 



A(m|,m 2 ) = 
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Fig. 2: The minimal fine-tuning as a function of the messenger scale M$ for tan (3 = 
10. The top black line is the total minimal fine-tuning as defined in equation (jlOj) 
which includes all the individual contributions. The individual contributions to 
the fine-tuning from fi 2 , rrijj , the gaugino masses M 2 , Mf and Mf, and the stop 
soft trilinear coupling A 2 are included. Moreover, the average fine-tuning of the 
stop soft masses and m \ R is included as in equation (jlip . 



The (top) black line shows the total FT as denned by (jTOj) . 

From the plot it is clear that the MFT increases as a function of the 
messenger scale Ms- This is expected since a higher messenger scale implies 
more RG running to the low scale so that small differences in high-scale 
input parameters are magnified. For Ms = Mgut, the total MFT is about 
22, i.e. 4.5%. (As an aside, for tan (3 = 30 and m A = 1000, the MFT for 
a Higgs mass of 114 GeV is about 11, i.e. 9%.) The largest contribution to 
the total minimal FT comes from Mf and A 2 which are both comparable for 
all values of Ms- The next most important contribution is that from Mf . 
The contributions from well as m~ L and m?^ are less important and 

increase only slightly as a function of Ms- The FT from m 2 Hu is very small 
for all messenger scales while the contribution from M\ is negligible for small 
and large Ms but larger for intermediate messenger scales. 

The large contribution from Mf is mainly because it has the largest (in 
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Fig. 3: The messenger scale values of M3, M2, Mi, A t and the average of the stop 
soft masses squared, m^, that give the minimal fine-tuning (MFT) as a function 
of the messenger scale Ms and for tan/3 = 10. The high-scale values of M2 
and At, and to a lesser extent Mi and m^, in the minimal fine-tuned region are 
roughly constant. The high-scale value of M3, however, decreases significantly as 
the messenger scale is increased. The reason for this is that the coefficient of M 2 
in the expression for m 2 z increases as a function of Ms, and thus the minimal 
fine-tuned region requires the value of M3 to decrease as Ms increases. 



magnitude) coefficient in the expression for m 2 z , at least for Ms > 10 10 GeV, 
see Figure [TJ The coefficients of the cross-terms A t M 3 , M 2 M 3 and M1M3 
are smaller (see Appendix [B]) , but together still contribute about 40% of the 
FT with respect to M| for Ms = Mqtjt- The reason that the cross-term 
contributions are so large is that the MFT values of A t , M2, and Mi are 
rather sizeable at the messenger scale when compared with M3 (at least for 
M s > 10 4 GeV). This is depicted in Figure H 

The FT of m| with respect to A 2 is also very large even though the 
coefficients of A 2 and the cross-terms A t M 3 , A t M 2 and A t M\ in the expression 
for m 2 z are rather small (for Ms = Mqut ; about 50% of the FT comes from 
the cross-terms). This is again because A t , M 2 and Mi are sizeable at Ms- 
The contribution to the FT from Mf is large for similar reasons. 

The FT with respect to fi 2 increases only slightly as a function of M5 
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-610 GeV 


305 GeV 


110 GeV 


475 GeV 



Table 1: Low-scale values for the stop soft trilinear coupling, the average of the 
left- and right-handed stop soft masses and the two physical stop masses. These 
low scale values give the minimal fine-tuning for arbitrary messenger scales. 



since the coefficient of /i 2 in the expression for m z does not vary much, and 
since the high-scale value of /i 2 increases only slightly as Ms is increased. The 
contribution from /x 2 is smaller than those from Mf , Mf and A 2 because the 
value of n is comparatively small and also because there are no cross-terms in 
the FT expression that involve /j, and other (large) soft parameters. Similar 

reasoning holds for the contributions from , m~ and m~ . 

to t L t R 

The low-energy spectrum that gives the MFT for a given messenger scale 
remains roughly unchanged as the messenger scale changes. The value of the 
stop soft trilinear coupling at the low scale is always about -610 GeV, with 
the two physical stop masses around 110 GeV and 475 GeV, respectively, 
see Table [1] and Figure HI These values of the stop-sector parameters are 
essentially determined by the constraint on the Higgs mass and from the 
minimization of A(m|, m 2 Hu (Ms)). The ratio X t /rrii is approximately -2, 

where X t = A t — /i cot /3, and m t - = \J\ijn 2 L + m i R )- The MFT is thus found 
for the natural maximal-mixing scenario which approximately maximizes the 
radiative corrections to the Higgs sector for a given set of parameters and for 
negative A t [3U HU H5l H6] . Small deviations of A t (and to a lesser extent 
and mf ) from its MFT value at the low scale lead to a very large increase in 
the FT, mainly from A(m|, m 2 Hu (Ms)). This can be seen from (I23p . which 
shows that the largest coefficients in the expression for m 2 Hu (M) in terms of 
low-scale parameters all involve powers of A t . Note that for generic points in 
the still allowed parameter space, A(m|,mjj ) would give one of the largest 
contribution to the FT. To minimize the FT it is thus best to minimize 
A(m|, m 2 Hu (Ms)) which essentially determines the values of the stop-sector 
parameters (see the discussion in Section f4 . 3 p . The other contributions to 
the FT are then not at their minimum, but they are much smaller and less 
sensitive to variations in the parameters. 

The low-scale values of the gaugino masses that give the MFT for a given 
messenger scale are shown in Figure HI While the value of M 2 that gives the 
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Fig. 4: The low-scale values of the gaugino masses Mi, Mi and Ms, the stop soft 
trilinear coupling A t and the average of the stop soft masses squared m 4 ~ that give 
the minimal fine-tuning (MFT) for the messenger scale Ms (with tan/? = 10). 
While the low-scale values of M2, At and that give the minimal fine-tuning are 
roughly the same for all Ms, the values of Mi and M3 decrease for larger Mg- 



MFT is roughly the same for all Mg, the values of Mi and M 3 decrease for 
larger Ms- Changing Mi away from its MFT value does affect the FT but 
not excessively so, while a change in M 3 has a larger effect. The /z-parameter 
is always found to be less than 150 GeV for the MFT region at any messenger 
scale. Choosing it to be closer to 100 GeV instead has a negligible impact on 
the FT, and allows a neutralino to be the lightest SM superpartner (LSP), 
instead of the lighter stop, which is found to be the LSP in the numerical 
minimization procedure. 

Negative A t may be expected to lead to less FT than positive A t because 
At has a strongly attractive infrared quasi-fixed point near [17J HH] 

A t ~ -M 3 . (12) 

(This relation is strictly valid only at the Pendleton-Ross quasi-fixed point 
for the top Yukawa [19], and neglecting SU(2)l and U(l)y gauge interac- 
tions.) Because of this it is most natural for A t and M 3 to have opposite 
sign and be comparable in magnitude at low scales due to renormalization 
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Fig. 5: The RG-evolution of At/M 3 for various low-scale boundary conditions 
A t (mz)/M 3 (m z ) = {-2.0, -1.5, 1.5, 2.0} and tan/3 = 10. The strongly at- 
tractive infrared quasi-fixed point near At/M% ~ — 1 is clearly visible. The gaugino 
masses have been set to their minimal fine-tuned values for the case Ms = Mgtjt> 
i.e. M 3 (m z ) ~ 335 GeV, M 2 (m z ) ~ 430 GeV, and M x (m z ) ^ 830 GeV. 



group evolution, see Figure [5j For positive A t and maximal-mixing in the 
stop-sector, A t would have to be an order of magnitude larger then M 3 at 
the messenger scale (see Figure |5]) which would lead to a much more FT 
parameter region. The MFT region here does not satisfy ([12]) exactly, but 
instead A t /M 3 ~ —1.8 at the low scale, for Ms = Mqut- In order to satisfy 
( |T2|) exactly, M 3 would have to be larger (assuming A t remains fixed). This 
would increase the size of the stop masses under RG evolution as can be seen 
from their /3-functions, see (J50l) and (JBTT) . which would lead to increased FT. 

The MTMSSM has negative soft squark squared masses at the messenger 
scale (see also [19] ) . This remains the case even if the messenger scale is very 
low and only on the order of a few TeV (for very low messenger scales, finite 
threshold corrections should really be included). Under RG-evolution the 
masses get driven positive very quickly within about a decade of running. It 
is the sizeable values of the gaugino masses that pull them up towards positive 
values. For smaller messenger scales the MFT region has a larger gluino mass, 
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Fig. 6: The RG-trajectories of the minimal fine-tuned region if the messenger scale 
is Ms=Mqut (tan/3 has been set to 10). At the scale mz, the parameter values 
are m~ t ~ 305 GeV, m k ~ 110 GeV, m h ~ 475 GeV, M 3 (m z ) ^ 335 GeV, and 
^j(m^) = 140 GeV. The minimal fine-tuned value is obtained for natural maximal- 
mixing, i.e. At — 2m,£. 



which drives the squark masses to positive values even faster while running 
towards the infrared. Equations (1221) and (124j) or (1251) in Appendix |A] show 
that negative squarks at the messenger scale lead to more stop-mixing at the 
low scale, as was pointed out in [19]. Figure [6] shows the RG-trajectories of 
the MFT region if the messenger scale is Ms=Mgtjt- 

The presence of tachyonic squarks at the messenger scale [501 EI] and/or 
very large A t [52| 153] may lead to dangerous color and / or charge breaking 
(CCB) minima. 

Very large A t may result in dangerous CCB minima around the EW scale. 
These CCB minima occur in the (£l,£r, H u ) plane [51] • The condition that 
the EW minimum is the global minimum may be estimated by going along 
the D-flat direction \ti\ = \tp\ = \H U \ and is given by [55] 

A 2 t +3fi 2 < 3(ml+ mfj. (13) 

Assuming instead that the EW minimum is only metastable but has a large 
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enough lifetime gives the weaker constraint [SS] 



A 2 t +^ 2 < 7.5(mf +m? ). 



(14) 



The MTMSSM easily satisfies the second condition, as well as satisfying the 
first condition. There are thus no dangerous CCB minima resulting from 
large A t . 

Tachyonic stops at the messenger scale may result in an unbounded from 
below potential along D-flat directions involving the stop fields, as well as 
first and/or second generation squark fields or slepton fields. Loop corrections 
give rise to an effective potential which is not unbounded from below, but 
they generically introduce a CCB minimum with a vacuum expectation value 
(VEV) on the order of the messenger scale. The MTMSSM may thus have 
CCB minima with a VEV around the EW scale if the messenger scale is 
low, or CCB minima with a VEV large compared to the EW scale if the 
messenger scale is high. Since the EW minimum is metastable and long-lived 
for mi > |M 3 [56J, it turns out that these CCB minima are not dangerous 
in the MTMSSM. Moreover, the MTMSSM does not determine the masses 
of the sleptons or first and second generation squarks since these do not play 
an important role in the FT. It is thus always possible to choose them in 
such a way to avoid CCB minima without changing the above FT results. 

Finally, it is interesting to note that there are several near degenerate 
parameter subspaces along which the FT does not change much. The first and 
second generation particles and their superpartners do not contribute much 
to the FT because in equation (j3j) they appear only with a small coefficient. 
The parameter Sy is also not very important for the same reason. A more 
interesting near degenerate subspace is that the FT is rather insensitive to 
changes in the difference of the two stop soft mass squared parameters at the 
low scale as long as their sum is kept fixed. This may be understood from 
the expression for m|, e.g. equation (J7|), in which only their sum appears 
(using the one- loop RG equations). However, even with only one-loop RG 
equations this degeneracy is not exact since small discrepancies appear in 
the FT measure from equations (1241 and ( |25l) . Moreover, the difference in 
the two stop soft mass squared parameters appears in the calculation of the 
physical stop masses and this affects the size of the Higgs mass, which is the 
most crucial low-energy constraint when calculating the FT. The FT only 



starts to change by an order one number when . / |m~ — m~ I ~ 300 GeV for 




M s = M GUT . 



17 



4.3 Analytic Motivation for Numerical Results 

The numerical results presented in section I4~2l may be motivated analytically. 
The discussion will for now assume Ms = Mqut, but generalizes to arbitrary 
Ms with a few caveats discussed below. 

In order to get a physical Higgs mass satisfying the experimental bound 
without generating large FT for the EWSB, it is natural to maximize the 
radiative corrections to m^. Due to the strongly attractive quasi- fixed point 
for At, this is achieved for negative A t near (natural) maximal mixing (at 
least for rrih not too large, see Section E}. 

The most important contribution to the FT comes from A(m|, m 2 Hu (Ms)) 
since it has the largest coefficients, see Appendix [B] Eliminating rh 2 H with 
the EWSB equation ([7]) and using the average stop soft mass squared m? = 
(m| + m| )/2 gives 

m%A(m%,m 2 H J = | - m% - 2.19 fi 2 + 1.36m? + 5.24 M 2 

-0.44M 2 + 0.46M 3 M 2 - 0.77 A t M 3 - 0A7A t M 2 
-0.01M 2 + 0.22 A 2 \. (15) 

It is possible to have cancelations among the various terms in this expression. 
A(m?z, M 2 (Ms)) also has large coefficients, but cancelations among its terms 
are impossible since A t is negative (see Appendix IB]) . 

Ignoring /a 2 , cancelation of the largest terms in equation ffTBI) . i.e. the 
gluino term and the average stop soft mass squared term, decreases the FT 
by setting m 2 H ~ m 2 H and leads to tachyonic squarks at the messenger scale 

m 

m\ ~ -3.9M 3 2 . (16) 
Next, the four terms on the second line of equation ffTSI) can cancel by taking 

M s K a96 ^ + 37i '. (17) 

Assuming M 2 — — A t , this simplifies to M 2 — 4.5M3. Furthermore, keeping 
only the most important terms, the natural maximal-mixing scenario implies 



~— ~ (0.32i t - 2.13M 3 - 0.27M 2 - 0.03Mi) 



0.66m? + 5. 15M 3 2 
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-1/2 

-OHM 2 + 0.02M 1 2 + 0.19A t M 3 + 0.04 A t M 2 - 0.05A 2 



;-4.80M 3 - 0.03M X ; 



1/2 

2.I6M3 2 + 0.02M 2 (18) 



which leads to M\ ~ I5M3, again assuming Mi — —A t . It is now possible to 
compute the ratio of the soft trilinear coupling with the gluino mass at the 
EWSB scale, 

A t 0.32i t - 2.13M 3 - 0.27M 2 - 0.03Mi , , 

— ~ ~ —1.8. 19) 

M 3 2.88M3 

These results agree well with the numerical results presented in section 14.21 
Note that a GUT scale model which predicts degenerate and negative 
squark and slepton soft masses at the GUT scale would need very large wino 
and bino masses in comparison to the gluino mass in order to drive the 
slepton soft masses to positive values under RG running to the EWSB scale 
[57] . This is due to the small coefficients of the bino and wino masses in the 
/3-functions of the slepton soft masses. It is interesting that the MFT region 
prefers the bino mass larger than the wino mass and, in turn, the wino mass 
larger than the gluino mass. 

Although this cancelation pattern holds to a good approximation for 
higher messenger scales, m~ does not exactly cancel Mf as the messenger 
scale decreases. For lower messenger scales, m~ becomes less tachyonic while 
Mf increases, allowing the stop masses to be driven positive faster under RG 
running to the EWSB scale. Moreover, the coefficient of Mf in the expres- 
sion for m 2 z ([5]) decreases significantly, as can be seen in Figure [TJ Therefore 
the cancelation pattern in A(m|,m 2 ^ ) discussed above does not hold since 
the m| contribution decreases while the Mf term gives a comparable con- 
tribution for all messenger scales (except for very small messenger scales). 
On the other hand, being a supersymmetric parameter, fi and its coefficient 
in equation (jSJ) does not change much for different messenger scales. Com- 
pared to Mf and m|, its contribution becomes important at lower messenger 
scales and a lower FT can be obtained by canceling the three contributions 
together. The other relations in the above cancelation pattern holds to a 
good approximation for lower messenger scales, although for M$ ^ 10 5 the 
cancelation pattern becomes more involved. 
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4.4 Summary of Phenomenological Implications 

The above analysis shows that the MTMSSM has small values for /i, the 
stop masses and the gluino mass. The gluino in the MTMSSM is around 
335 GeV for Ms = Mgut, but heavier for lower Ms- There is large mixing 
in the stop-sector which introduces a significant splitting between the two 
physical stop masses. They have masses of around 115 GeV and 475 GeV 
respectively, see Table [TJ Thus the MTMSSM may have a stop as the LSP. 
However, as mentioned before, /i can be chosen to be small enough so that a 
neutralino is the LSP without affecting FT by much. 

At the Large Hadron Collider, gluino pair-production in the MTMSSM 
is thus rather large and comparable to top quark pair-production. The pro- 
duction of Mi is also of the same order. 

The gluinos are Majorana particles, and can decay into the lightest stop 
via gg — > tttiti producing same-sign top quarks 50% of the time. The top 
quarks each decay into Wb, and the events with two same-sign top quarks will 
contain two same-sign leptons if the W decays leptonically. If a neutralino 
and a chargino are lighter than the stop, the decay t\ — > xfb is possible, with 
Xi further decaying into a neutralino and soft jets or leptons. The events 
thus also contain missing energy and a number of 6-jets, some of which are 
soft if the ti — xt mass splitting is small. 

If ti is the LSP a number of further interesting signatures are possible, 
see [SH]. The lighter stop can either be pair-produced directly or from gluino 
decays. Even though it is the lightest SM superpartner, it may decay into a 
lighter goldstino G via the flavor- violating decay t\ — > cG or via the three- 
body decay ti —+ bWG. The decay rate depends on the messenger scale, 
with lower messenger scales leading to larger decay rates. For reasonable 
messenger scales, its decay length easily exceeds the hadronization length 
scale, and the stop in general hadronizes before it decays [58J. For messenger 
scales less than a few hundred TeV, the decay length is small enough so 
that the decay products seem to originate from the interaction region. The 
three-body decay leads to a similar signature as the top decay but can be 
distinguished from it, see [5H]. For larger messenger scales, ti decays inside 
a hadronized mesino or sbaryon and a variety of interesting signatures are 
possible (58], including mesino-anti- mesino oscillations |60j . 

Another interesting possibility is the direct pair-production of the heavier 
stop t 2 . Since the two physical stop masses are split by a large amount, 
the decay mode t 2 — ► ti + Z is kinematically allowed and has a sizeable 
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branching ratio [M] • The resulting signature depends on the t\ decay channel 
as discussed above. For \t an d Xo lighter than t\, the authors of [61J propose 
to look for the inclusive signature Z(l + , l~)bbpxX , where the two leptons l + 
and l~ have an invariant mass equal to the Z-mass. Detecting this signature 
would give evidence for the maximal- mixing scenario but requires a large 
integrated luminosity (at least (9(100 fb -1 )) [5T]. Since the mass difference 
between t\ and the LSP is small in the MTMSSM this signature will be very 
hard to see since the jet from the decay t\ —* xtb is soft which makes it more 
difficult to separate the signal from the SM background [61J. 

An alternative way to measure the parameters in the stop-sector is to use 
the Higgs boson as a probe [02] • A measurement of the Higgs mass and its 
production rate in the gluon fusion channel allows the average of the two 
stop soft masses as well as the stop mixing to be determined in many regions 
of the still allowed MSSM parameter space, and especially in regions where 
the FT is small [62]. 

4.5 Fine-Tuning with Respect to Other Parameters 

This subsection briefly discusses other parameters that may in principle con- 
tribute to the FT. 

If the goal is to find the MFT region of a model and make a prediction of 
what parameter region is preferred for the model from a FT point of view, 
there is no reason to include the FT of experimentally known parameters 
such as gy, #2 ,<73> or X t . Taking into account the known parameters in the 
minimization procedure would most likely lead to other MFT values for all 
parameters, including MFT values for the known parameters which would in 
all likelihood not match the experimental values. 

If the goal, however, is to find the FT of a given model, one should in 
principle include contributions from experimentally known parameters. For 
example, FT with respect to A t , A(m|, X t (Ms)), may give a large contribu- 
tion to the total FT due to the large top mass. Indeed, with the MFT values 
for M s = M GU t, A(m|, At(Mftf GUT )) — 8. This, however, increases the total 
FT only by a small amount from 22.1 to 23.5. 

What about FT with respect to m\ 2 and tan/5? These parameters are 
unknown and in principle they should be included in the minimization pro- 
cedure. With the help of equation ([3]) and symmetries, it is however easy to 
see that A(m|, m\ 2 (M s )) = 0. Indeed m\ 2 does not appear directly in the 
expression for m 2 z . Furthermore it breaks a [/(I)pq- and a i7(l ^-symmetry 
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and consequently does not feed back into any other /3-functions since no other 
parameter breaks both symmetries. Thus m\ 2 cannot appear in equation (jSJ) 
and is therefore completely free, which allows to be chosen accordingly 
as discussed in Section [2j 

The FT of tan/3 has not been taken into account in the minimization 
procedure since an explicit expression for m| can only be obtained assuming 
a specific value for tan /3, because A t depends on tan (3 through m t . Moreover, 
since tan (3 is then a free parameter the approximation leading to equation (111 
may not be valid anymore and m 2 H should be reintroduced. Contributions 
from bottom/sbottom and tau/stau sectors should also be included if tan/3 
becomes large. 

5 Minimal Fine- Tuning as a Function of the 
Higgs Mass 

The Higgs mass is the most important low-energy constraint that deter- 
mines the amount of minimal fine-tuning (MFT). It is therefore interesting 
to look at how the MFT is affected when the lower bound on is changed. 
Figure [7] shows a plot of the MFT as a function of the lower bound on m^, 
where the calculation of rrih is the same one used in the FT minimization 
described in Section 14.11 and only includes the one-loop corrections from the 
top-stop sector (with = 250 GeV, tan/3 = 10, m t = 170.9 GeV, and 
Ms = Mqut)- The Higgs mass calculated with the one- loop corrections will 
be denoted by mjf. The region of MFT always saturates the bound on m\ e 
and has negative A t . The minimal FT is about 1% for mjf ~ 132 GeV. 

There are, however, other important one-loop and two-loop corrections 
that can significantly affect nth, and these need to be included in order to 
get a more accurate idea of how the MFT chang function of the 

lower bound on m^. With these additional corrections, is not anymore 
a symmetric function of the stop- mixing parameter X t = A t — //cot (3 ~ A t , 
where the latter approximation is good for sizeable tan (3. It can be up to 5 
GeV larger for X t = +2m 4 - than for X t = — 2m 4 -, the difference arising from 
non- logarithmic two-loop contributions to rrif L , see [63| [Ml ES]- Moreover, 
large chargino masses, i.e. large values of M 2 and /i, can give important 
negative contributions to [66]. These corrections are also not included in 
m\ e . Two-loop corrections that allow the gluino mass to affect m/, can also 
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Fig. 7: The minimal fine-tuning as a function of the lower bound on the Higgs 
mass to/j, where the calculation of only includes the one-loop corrections from 
the top-stop sector (tan/3 = 10, ttia = 250 GeV, m t = 170.9 GeV). 



be important but are smaller in general - this will be ignored in the following 
discussion since the impact on the results presented below is negligible. 

The MFT spectrum that was found with the minimization procedure may 
be used to calculate with FeynHiggs. The FeynHiggs estimate for 
will be denoted by m™. The result is the solid black line in Figure [BJ This 
MFT spectrum characteristically has large chargino masses and a negative 
value for A t near the "natural" maximal mixing scenario. 

Comparing the solid black line in Figure[8]with the curve in Figure [7] shows 
the well-known fact that the higher-order corrections to rrih are extremely 
important. There are two additional very striking features. First of all, as 
m™ increases and approaches 120 GeV, the FT increases enormously. Any 
further small increase in the Higgs mass results in an exponentially large 
increase in the FT. The reason is that as m™ approaches 120 GeV here, it 
only grows logarithmically as a function of the stop masses. The stop masses 
therefore become exponentially large and thus increase the FT exponentially 
(see also [31]). 

The second striking feature of this curve is that the value of the Higgs 
mass at which the FT starts to increase exponentially is rather low (the MFT 
is already 1% for m™ ~ 119 GeV). This value of rrih may be increased by just 
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under 2 GeV by choosing larger tan /3 and rriA (recall that throughout this 
discussion tan/5 = 10 and = 250 GeV). Note that the latest Tevatron 
top mass value (m t = 170.9 GeV) has been used in the calculation, and a 
slightly different value can also change nih by a few GeV. 

An obvious question is whether the MFT region is significantly different 
if m™ were used in the minimization procedure instead of m\ e (the former 
is too complicated to be used). For MSSM spectra that give small rrih this 
is certainly not the case, since there is not a very large discrepancy between 
the two Higgs mass estimates mjf and m™. The difference between the two 
Higgs mass estimates becomes significant, however, for MSSM spectra that 
give a large rrih, and the approximation m™ can be substantially smaller than 
m h e . Also, as mentioned above, m™ can be substantially larger for positive 
A t (near maximal mixing) than for negative A t (near "natural" maximal 
mixing), and increases as the chargino masses decrease. On the other hand, 
mjf remains unaffected by the sign of A t and the size of the chargino masses. 
It is thus possible that the MFT region does not coincide with the region 
obtained in the above minimization procedure as the lower bound on 
increases. This is indeed the will now be discussed. 

The FT may be minimized with the constraint that the chargino masses 
are small. Since the effect of varying fi and M2 on the FT are noticeable 
but not substantial, the resulting spectrum will be characterized by gluino 
and stop masses that are only slightly larger than those obtained in the MFT 
region discussed in this paper. The value of A t is still negative. This spectrum 
may be used to calculate m F h E . The result is shown by the dash-dot green 
curve in Figure [HJ For m™ not too large, the solid black curve lies below the 
dash-dot green curve because the MFT region has large values of M2, see 
Section HI As m™ increases further, however, the FT becomes exponentially 
large since the stop masses become exponentially large. Smaller chargino 
masses lead to larger values of m F h E , and the two curves show that for rrih just 
below 120 GeV, a smaller FT may be obtained by decreasing the size of M2. 
This behavior cannot be captured by m h which is unaffected by a change 
in the chargino masses. Note that the transition between the two regions 
described by the two curves is smooth, and that it occurs when the MFT is 
already more than 1%. 

Next, the FT may be minimized with the constraint that A t is positive and 
near maximal mixing. The resulting low-energy spectrum is characterized 
by small chargino and gluino masses. This spectrum may then be used to 
calculate and the MFT as a function of this value of m™ is displayed 
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Fig. 8: The minimal fine-tuning as a function of the lower bound on the Higgs 
mass m/j calculated with FeynHiggs 2.6.0 (tan (3 = 10, rriA = 250 GeV, mt = 170.9 
GeV). Throughout this paper the fine-tuning is minimized subject to a constraint 
on rrih, where rrih is estimated with a one- loop formula as described in Section T4.ll 
The different lines arise from different assumptions made about At, or fj, and M2, 
when minimizing the fine-tuning. These different assumptions give rise to different 
low-energy spectra that present the least fine-tuned parameter choices satisfying 
these assumptions. These low-energy spectra may then be used in FeynHiggs 
to calculate rrih- Although M2, fi and the sign of At do not affect the one-loop 
estimate of rrih which only contains the dominant corrections, they do affect the 
FeynHiggs estimate of rrih- For the solid black line no constraint was set on At, 
and \x and M2 were only required to be above 100 GeV. It is the same line as in 
Figure but with rrih estimated by FeynHiggs instead of the one-loop formula. 
The dashed blue line assumes At is positive and near maximal mixing, also with 
M2 and \i only required to be above 100 GeV. The dash-dot green curve makes 
no assumption about A t but sets fi = 100 GeV and M2 = 100 GeV. The dotted 
red line assumes At = 0, and again only requires fi and M2 to be larger than 100 
GeV. Further details and explanations are given in the text. 
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by the dashed blue line in Figure [SJ Comparing the solid black line or dash- 
dot green line with the dashed blue line, it is clear that for small m™ the 
MFT region has negative values of A t . Even though negative A t might be 
expected to always give less FT than positive A t due to the IR quasi-fixed 
point, the increase in m™ by several GeV by making A t positive is substantial, 
and as m™ approaches about 123 GeV, the two curves cross. Thus, there 
is a transition from A t ~ — 2m^ to A t ~ +2m^ of the minimal fine-tuned 
region as m™ increases. This behavior is again not captured by m\ e which is 
independent of the sign of A t . The transition occurs when the minimal FT 
is already quite large (about 0.2%). 

This transition from negative to positive A t is not smooth, in the sense 
that the first derivative of the curve at the transition point is not continuous^]. 
To show this, the FT may be minimized with the constraint A t — 0. The 
resulting low-energy spectrum may then again be used to calculate and 
the result is shown by the dotted red line in Figure [H The value of m™ 
for vanishing stop-mixing, A t = 0, is much lower than for the two maximal 
mixing scenarios, A t ~ ±2m 4 ~, and it is clear that the MFT region does not 
interpolate smoothly between them as a function of A t . 

The main point of the analysis in this section is that although the MSSM 
is already fine-tuned at least at about the 5% level (if the messenger scale 
equals the GUT scale), there is not much room left for the Higgs mass to 
increase before the FT becomes exponentially worse. 

Note that for a lower messenger scale the Higgs mass can have a slightly 
larger value before the MFT begins to increase enormously. For example, 
for M s = 200 TeV, the MFT is 1% for m h ~ 123 GeV. So even for a lower 
messenger scale the Higgs mass cannot be that much beyond 120 GeV before 
the MFT increases dramatically. 

6 Conclusions 

This paper presented the minimally tuned Minimal Supersymmetric Stan- 
dard Model (MTMSSM). The MSSM parameter region that has the minimal 
model-independent fine-tuning (FT) of EWSB was found. Model - indepen- 
dent means that no relations were assumed between the soft SUSY breaking 
parameters at the scale at which they are generated, here referred to as the 
messenger scale. Instead, all of the important parameters were allowed to be 

2 One may perhaps refer to this as the first order phase transition of fine-tuning. 
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independent and free at the messenger scale, and were taken to contribute 
to the total FT of the EWSB scale. The messenger scale itself was varied 
between 2 TeV and Mqut and the effect of this on the minimal FT was 
presented. 

The most important parameters that contribute to the tuning are |//| 2 , 
m 2 Hu , the gaugino masses Mi, M 2 and M 3 , the stop soft masses and , 
and the stop soft trilinear coupling A t . The MSSM spectra which lead to the 
minimal model-independent FT were found by numerically minimizing the 
FT expression subject to constraints on the Higgs, stop, and gaugino masses 
(the Higgs mass was found to always be the most important low-energy con- 
straint). The high-energy spectra are characterized by tachyonic stop soft 
masses, even for messenger scales as low as 2 TeV (but note that threshold 
effects in the RG- running were neglected throughout). The potential exis- 
tence of charge and/or color breaking minima turns out not to be a problem. 
The gluino mass, M3, is much smaller than the wino mass, M2, and M2 in 
turn is much smaller than the bino mass M\. The low-scale spectra are char- 
acterized by negative A t near the maximal mixing scenario that maximizes 
the Higgs mass. The large stop mixing leads to a large splitting between the 
two stop mass eigenstates. Interesting phenomenological signatures include 
the possibility of a stop LSP. 

The minimal FT was also found as a function of the lower bound on 
the Higgs mass (with the messenger scale set to M GU t)- Although in the 
numerical minimization procedure the dominant one-loop expression for rrih 
was used as a constraint, the resulting least fine-tuned spectra were used to 
calculate rrih more accurately with FeynHiggs. A plot of the minimal FT as a 
function of rrih was presented. There are several striking features of this plot. 
For rrih larger than about 120 GeV the FT increases very rapidly. This value 
of rrih is rather low, perhaps surprisingly so. It is only slightly dependent on 
the parameters in the Higgs sector. Near it, the value of A t in the least FT 
region also makes a sudden transition from lying near — 2mj to lying near 
+2m£, where m 4 ~ is the average of the two stop soft masses. The upshot of 
this particular analysis is that although the MSSM is already fine-tuned at 
least at about the 5% level (if the messenger scale equals the GUT scale), 
there is not much room left for the Higgs mass to increase before the FT 
becomes exponentially worse. 
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A Semi-numerical Solutions of the MSSM 
One-Loop RG-Equations 

This appendix reviews the procedure for solving the MSSM one-loop RG 
equations semi-numerically [351 [36] • The low scale Mq is set to be mz, and 
the high (messenger) scale Ms is taken to lie anywhere between mz and 
Mgut- Threshold corrections are neglected when solving the RG-equations. 

The main goal is to obtain an expression for m 2 z in terms of high-scale 
input parameters as in equation (J5J). Assuming that tan/3 is not too small, 
this requires solving |/i(m^)| 2 and m 2 Hu (mz) in terms of high-scale parameters 
(for moderate values of tan/3, mfr may be neglected, see equation (J4])). The 
fine-tuning may then be calculated and naturally expressed in terms of high- 
scale parameters as in equation fllOp . However, in order to minimize the 
fine-tuning taking into account low-scale constraints on the Higgs, stop, and 
gaugino masses, it is more appropriate to rewrite the fine-tuning expression 
in terms of low scale parameters. This requires that \i as well as all the soft 
supersymmetry breaking parameters appearing in equation (flOl) be written 
in terms of low scale parameters. 

In solving the RG-equations, only the contributions from the third gener- 
ation particles will be included, since the third generation Yukawa couplings 
are much larger than those from the first and second generations. More- 
over, the contributions from the bottom/sbottom and tau/stau sectors are 
neglected as tan /3 is taken to be not too large. 

The high-scale parameters may in general be written in terms of low 
scale-parameters as 

m 2 (M s ) = Cijfc(tan/3, M , M s ) m,(M ) m k (M ). (20) 
For example, for Ms = Mqut ; the expressions for the most important high- 
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scale parameters written in terms of low-scale parameters are 

Mi = diMi {d 1 ,d 2 ,d 3 } = {2.42,1.22,0.35} (21) 

A t = 3.15 A t + 2.33 M 3 + 1.03 M 2 + 0.26 M 1 (22) 
m 2 Hu = 2.07 m 2 Hu + 1.07 + 1.07 mf fl + 0.19 Af 3 2 - 0.98 M 2 

-0.31M 1 2 + 3.38 A 2 + 3.69 A t M 3 + 1.19A t M 2 + 0.24 A f Mi 
+ 0.76 M 3 M 2 + 0.15 M 3 M x + 0.05 M 2 Mi + 0.06 S y (23) 
m} L = 0.36 m 2 ^ + 1.36 + 0.36 m\ R - 0.72 M\ - 0.81 M 2 2 



- 0.06 Mf + 1.13 + 1.23 A t M 3 + 0.40 A t M 2 + 0.08 A t M x 

+ 0.25 M 3 M 2 + 0.05 M 3 M 1 + 0.02 M 2 Mi + 0.02 S Y (24) 
m, 2 = 0.72m 2 + 0.72m 2 + 1.72m 2 — 0.65 M? — 0.18 M| 

- 0.46 M 2 + 2.26 A 2 + 2.46 A t M 3 + 0.80 A t M 2 + 0.16 A t Mi 

+ 0.50 M 3 M 2 + 0.10M 3 Mi+ 0.04M 2 Mi- 0.09 S Y (25) 
fi = 0.95 p. (26) 

Similar type of expressions hold for low-scale parameters as a function of 
high-scale parameters. The gauge couplings g a , a G {1,2,3}, and the top 
Yukawa coupling At are fixed at the low scale by their experimental values 
[40] . Section [A.ll gives the solution of their RG-equations. 

The MSSM one-loop /3-functions that need to be solved come in three 
different functional forms [67]. The RG-equations of the gaugino masses M a , 
the supersymmetric Higgsino mass /i, and Sy are of the form 

din ■ 

"~j7~ = /i(At, g a ) rrii, m^ e {M a , /i, S Y }, (27) 
where t = \n(Ms/M ). Their solution is given by 

rrii(t) = mi(0) exp / dt' fi(X t ,g a ). (28) 

JO 



The stop soft trilinear coupling has the functional form 

a{X t )A t + b(g a ,M a ). (29) 



cL4t 



r/f 

The solution of this equation is more involved due to the presence of both 
homogeneous and inhomogeneous terms, and requires the solution for the 
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gaugino masses (1251) . It may be written as (see Section 1X73]) 

A t (t) = e f dt ' aiXt) A t (0) + ef dt ' aiXt) [ dt'e-f dt " a{Xt) b{g a ,M a ). (30) 



Finally, the RG-equations of the up-type Higgs soft mass and the stop 
soft masses form a system of coupled inhomogeneous differential equations, 

^ L = ^2,Uij{\t)m) + Vi(g aj M aj S Y ,At), m- G {m| u , mf^, mf fl }. (31) 



This may be solved (see Section IA.4j) using the solutions for the gaugino 



masses and Sy fl28l) as well as the solution for the stop soft trilinear coupling 

(EQD, 

m?(t) = ^ dt ' uiXt) m 2 {0) +ef dt '^ dt'e-f dt ''^v(g a ,M a ,S Y ,A t ) 



(32) 



A.l Gauge and Yukawa Couplings 



The one-loop /^-functions for the gauge and top Yukawa couplings in the 
MSSMare 

Svr 2 /^ = b a gi, {b Y , b 2 , b 3 } = {11, 1,-3} (33) 

167r 2 /5 At = A, ( 6 X 2 - f g\ - 3 g\ - H g ^ . (34) 

Their solutions are 

gl(t) = gl(o)C\t) (35) 

X 2 t (t) = X 2 t (0)E(t;n )G(t;n y\ (36) 

where n = ^, ^-J = (^§, 3, — y), and for future convenience the func- 
tions 

Ut) = l-^&W (37) 

3 

E(t;n) = H&>-(t) (38) 

a=l 
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F(t;n) = [ dt'E{t';n) (39) 
Jo 

G{t-n) = 1- A A 2( )F(t;n) (40) 



have been introduced. The solution (1361) is analytic if #2 and gy are set 
to zero [68l [69] , whereas non-zero values of #2 and gy require a numerical 
integration. 

A. 2 Gaugino Masses, /i-term and Sy 

The RG-equations for the gaugino masses, /i and Sy are 

A*. = ^T&S ( 41 ) 



16tt 2 ^ = /x (3X 2 -3g 2 2 -g 2 Y ) (42) 



The general solution is of the form (1281) . and may be written as 

M a (t) = M a (0)C\t) (44) 
fi(t) = K0)G{t;ti )-^l(t)t?(t) (45) 

Sy(t) = Sy(0)^\t) (46) 

with the notation of Section IA.ll The solutions for the gaugino masses and 
Sy are analytic while \x must be solved numerically unless the contributions 
from g 2 and gy are neglected. 

A. 3 Stop Soft Trilinear Coupling 

The /3-function of the stop soft trilinear coupling is 

8ir 2 (3 At = (d A 2 A t - y g 2 , M 3 - 3 g\ M 2 - ^ g 2 Y M^j . (47) 
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Using the solutions for the gaugino masses (I44|) . this equation may be inte- 
grated and written as 



G(t;n ) 



a=l 



(48) 

where (e°)/3 = are the usual unit vectors. If g 2 and gy are zero, the 
solution does not require a numerical integration. 



A. 4 Up-type Higgs Soft Mass and Stop Soft Masses 



The /^-functions of ml, , m~ and m~ are 

r H u 1 t L t R 

= 3Af 



87T 2 ^ m 2 



m 2 ^ + m~ r + m 2 n + |v4 t | 2 



-3g%\M 2 \ 2 - tilM^ 2 - ~&S Y 



(49) 



A 



m 2 M + m~ + m 2 + |A| 2 



2A 2 



16 

y 

2,2,2 



£3 |M 3 | 2 - 3 g 2 2 \M 2 \ 2 - l -g 2 Y |Mx| 2 - ^ g 2 Y S Y (50) 



1 

6 
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tie 



IAI 2 



f ^|M 3 | 2 -y^|Mi| 2 -^^y. (51) 

They form a system of coupled inhomogeneous differential equations. Note 
that At appears quadratically in these /5-functions which gives cross-terms 
between M a (0) and A t (0) (see equation (jUJ)). The equations can be solved 
as in f[3"2"j) but it is possible to simplify the analysis by the change of variables 



X 
Y 
Z 



WItt — rrir — Vtlr 



m 2 rr — 3m~ 

tlu t L 
2,2,2 

m M + rrir + m,~ 

LIU 



In terms of the new variables, the /3-functions are 
8ti 2 P x = ^ 3 2 |M 3 | 2 + ^ 2 ,|M 1 | 2 + ^ 2 S Y 



(52) 
(53) 
(54) 



(55) 
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16g 2 3 \M 3 \ 2 + 6g 2 2 \M 2 \ 2 -Ig*,^ 2 
6\ 2 t Z + 6\ 2 t \At\ 2 - 32 " 2 



g z 3 \M 3 \ z -6g z 2 \M 2 \ 
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(56) 

g*\M x \ 2 .{$7) 



In this form, fix and /3y are easily integrated since they have no homogeneous 
term (which is due to the fact that the corresponding matrix ity in ( 1311) has 
rank = 1) 



1 6 

X(t) = X(0)--M 2 (0) fe- 2 (t)-l) 



(58) 



+4 M i (°) (^ 2 w - !) + n s ho) (er x (0 - 1) 



y(t) 



m-gM^O) fe- 2 (*)-l) 



(59) 



1 



+3M|(0) (e 2 " 2 (t) - 1) - 33 (o) (er 2 (i) - 1) • 



The equation for Z requires a numerical integration (even if g 2 and gy are 
zero) 



Z(t) 



G(t; n ) 



m - E 



a=X 



n ) Q ^^ (<?(f ; n ) - £ (t) G(t; n - 2? 



+^\mj*dt'E(t'-,fi )\Mt')\ 2 
The solutions for mi, , m 2 and m~ in terms of X, Y and Z are then 



(60) 



X(f) + Z(i) 
X(t) - 2Y(t) + Z(t) 
±(-2X(t) + y(t) + Z(t)). 



(61) 
(62) 
(63) 



B Fine-tuning Components 



This appendix lists for completeness the expressions for the fine-tuning of 
m| with respect to Mf, M|, Mf, /i 2 , A 2 , m 2 Hu , ni~ and m| . The fine- 
tuning components as a function of high-scale parameters are easily found 
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from the fine-tuning measure, equation ([Hlh with the observable m 2 z written 
as in equation For M s = M GU t, the fine-tuning components are 



m|A(m|, Ml) ~ 5.24M 2 + 0.23M 3 M 2 + 0.03M 3 Mi - 0.38i t M 3 (64) 

m|A(m|, M 2 ) ~ -0.44M 2 + 0.23M 3 M 2 + 0.01M 2 Mi - 0.08itM 2 (65) 

m%A{m%, Ml) ~ -0.01M 2 + 0.03M 3 Mi + 0.01M 2 Mi - 0.0li t Mi (66) 

m|A(m|,/i 2 ) ~ -2.19/a 2 (67) 



m 



A(m|, A 2 ) ~ 0.22i 2 - 0.38i 4 M 3 - 0.08i t M 2 - 0.0L4 t Mi (68) 

m 2 z A(m z ,m 2 Hu ) ~ -1.32m^ (69) 
~ -m|- 2.19/t 2 + 1.36 m? + 5.24 M 3 2 

-0.44 if 2 + 0.46 M 3 M 2 - 0.77 A t M 3 - 0.17 A t M 2 
-0.01M 2 + 0.22 i 2 

m|A(m|,mfJ ~ 0.68m| (70) 

m|A(m|,m 2 fl ) ~ 0.68m 2 fl . (71) 



Here it is understood that the absolute value of the right-hand sides of each of 
these equations is meant to be taken. The EWSB relation, equation ([7j), was 
used to eliminate rh 2 H . It is natural to eliminate rh 2 H instead of fi 2 or any 
other soft supersymmetry breaking parameters since (x 2 is supersymmetric 
while the other soft supersymmetry breaking parameters are not involved in 
the EWSB equation at the EW scale. With the help of equations (!2T|) - (!26|) . it 
is now straightforward to rewrite the FT expression (fTUj) in terms of low-scale 
parameters. 
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